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A QCD inspired relativistic effective Hamiltonian model for the bound states of mesons has 
been constructed, which integrates the advantages of several QCD effective Hamiltonian models. 
Based on light-front QCD effective Hamiltonian model, the squared invariant mass operator of 
meson is used as the effective Hamiltonian. The model has been improved significantly in four 
major aspects: i) it is proved that in center of mass frame and in internal coordinate Hilbert 
subspace, the total angular momentum J of meson is conserved and the mass eigen equation can be 
expressed in total angular momentum representation and in terms of a set of coupled radial eigen 
equations for each J; ii) Based on lattice QCD results, a relativistic confining potential is introduced 
into the effective interaction and the excited states of mesons can be well described; iii) an SU(3) 
flavor mixing interaction is introduced phenomenologically to describe the flavor mixing mesons 
and the mass eigen equations contain the coupling among different flavor components; iv) the mass 
eigen equations are of relativistic covariance and the coupled radial mass eigen equations take full 
account of L — S coupling and tensor interactions. The model has been applied to describe the 
whole meson spectra of about 265 mesons with available data, and the mass eigen equations have 
been solved nonperturbatively and numerically. The agreement of the calculated masses, squared 
radii, and decay constants with data is quite well. For the mesons whose mass data have large 
experimental uncertainty, the model produces certain mass values for test. For some mesons whose 
total angular momenta and parity are not assigned experimentally, the model gives a prediction 
of the spectroscopic configuration 2 Lj. The connection between our model and the recent low 
energy QCD issues-the infrared conformal scaling invariance and holographic QCD hadron models 
is discussed. 

PACS numbers: 12.38.Le, ll.10.Ef, 14.40.-n 



I. INTRODUCTION 



To study hadronic properties at low energy scales, 
nonperturbative effects must be taken into account [l]. 
To describe mesons and baryons, there are several 
main approaches: coupled Bethe-Salpeter(BS) and 
Dyson-Schwinger(DS) equation approach, relativistic 
constituent quark model based on Bethe-Salpeter equa- 
tion(BSE), relativistic string Hamiltonian approach, and 
holographic light-front QCD approach. In the coupled 
Bethe-Salpeter and Dyson-Schwinger equation approach 
by P. Maris, P. Tandy, L. Kaptari et al.[2J|. the dressed 
quark propagators are assumed to have time like com- 
plex mass poles where the absence of real mass poles 
simulates quark confinement; the BS kernel is approxi- 
mated by ladder rainbow truncation with two-parameter 
infrared structure. The approach contains four parame- 
ters in u-d-s quark sector and is consistent with quark and 
gluon confinement. Besides, it has the feature of preserv- 
ing the relevant Ward identity and generating Dynami- 
cal chiral symmetry breaking. The vector mesons p, 0, 
and K* are studied in detail, the calculated masses of 
p, <^>,and K* mesons and decay constants f p , fy, and fa* 
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are within 5% and 10% of the data respectively. More- 
over, the ground-state spectra of light-quark mesons are 
also studied and a good description of flavor-octet pseu- 
doscalar, vector, and axial- vector meson spectrum is ob- 
tained. The applicable domain of ladder truncation and 
the relative importance of various components of the two- 
body BS amplitude are also explored. However heavy 
quark mesons are not investigated and the number of 
mesons treated are not too many. R. Alkofer, P. Watson, 
and H. Weigel Q follow the same approach, scalar and 
pseudoscalar, vector and axial vector mesons are studied. 
A similar approach is pursued by P. Jain and Munczekjj], 
about 50 mesons are investigated and the results are in 
good agreement with experiments. But heavy quarks are 
analyzed by non-relativistic dynamics. It should be noted 
that in contrary to Hamiltonian dynamics which works 
with wave functions that are not manifestly covariant 
quantities, the above BSE/DSE approaches emphasize 
the relativistic covariant aspect of the formalism and in- 
variant quantities are studied. 

The constituent quark model(CQM) works surpris- 
ingly well for most of the observed hadronic states [9, Q ■ 
However, as a phenomenological theory, there are still 
some problems and puzzles that need to be clarified and 
understood[7|. One of the most important problems is 
relativistic effect. To solve the relativistic covariant prob- 
lem of CQM, the relativistic constituent quark model 
based on Bethe-Salpeter equation was proposed by B. 



Metsch et al. in Bonn Group Q. In this approach, 
the meson and baryon Hamiltonians are extracted from 
Bethe-Salpeter equation and the relativistic covariant 
constituent quark models for mesons and baryons are 
constructed. Based on Dirac structure of the two-body 
effective interactions, two types of models( A and B) are 
constructed. This approach addresses hadron mass spec- 
tra from ground state to 3GeV, light-flavor mesons, scalar 
excitations, linear Regge trajectory, pseudoscalar mix- 
ing, and parity doublet (for baryons). In this approach, 
the Dyson-Schwinger equation(DSE) is approximated by 
parametrization of infrared effective gluon propagator, 
the interaction kernel of BSE is given by single gluon ex- 
change(OGE) and the confinement is parameterized by 
a string-like potential ( having two versions defined by 
Dirac structures A and B ). The instanton- induced spin- 
flavor dependent interaction is also included in the BSE 
kernel. The mass spectra up to 3 GeV, electroweak and 
strong-decay properties are calculated with 7 to 9 pa- 
rameters. About 60 scalar and pseudo-scalar, vector and 
axial vector, and some tensor mesons with J=0,l,2 are 
calculated by models A and B, and compared to Godfrey- 
Isgur's calculation and experimental data ( the deviation 
seems large but the errors are not indicated ). Due to 
the Dirac structure of the effective interactions, spin-spin 
and spin-orbital interactions are included. Besides, heavy 
mesons are not treated. 

The relativistic string Hamiltonian approach was pro- 
posed by A.M. Badalian et al.Q. The merit of this 
approach is that the quark-anti-quark interaction and 
the confinement are generated by the relativistic string 
( through Nambu-Goto action for QCD vacuum fluctua- 
tion) which leads to a large reduction of the number of 
model parameters. After quantizing the action by path 
integral, they construct a Hamiltonian with a linear con- 
fining potential and hyperfine quark-anti-quark interac- 
tions. Using only one parameter of string tension, they 
study the systematic property of orbital excitations and 
rotation of mesons. The linear Regge trajectory rela- 
tion between squared mass and orbital angular momen- 
tum is produced nicely and in agreement with the data 
for about 40 mesons. The relativistic string Hamilto- 
nian approach is spin-independent. In the lowest order, 
this approach doesn't contain spin-spin, spin-orbital, and 
tensor interactions, thus it can produce the spin averaged 
mass spectra for mesons. However, to include the higher 
order effects by perturbation method, the hyperfine spin- 
dependent interactions could be obtained. 

The holographic light-front QCD approach by S. J. 
Brodsky and G. F. de Teramond et al. [1(| is based 
on light-front QCD and AdS/CFT correspondence. The 
AdS/CFT correspondence between string theory in AdS 
space and conformal field theories in physical space-time 
leads to an analytic, semi-classical model for strongly- 
coupled QCD, which has scale invariance and dimen- 
sional counting at short distances and color confinement 
at large distances. This correspondence also provides 
AdS/CFT or holographic QCD predictions for the ana- 



lytic form of the frame-independent light- front wave func- 
tions (LFWFs) and masses of mesons and baryons. Re- 
cently, Brodsky et aZ.[lO| have found that the transverse 
separation of quarks within hadron is related to holo- 
graphic coordinate (the fifth dimensional z-coordinatc) 
in AdS/CFT correspondence, the mass eigen equation of 
meson in light-front effective Hamiltonian approach cor- 
responds to the equation of motion for the holographic 
field of effective gravity field of super string in AdS 
space at low energy limit. Recently, they have modi- 
fied the gravitation background by using a positive-sign 
dilaton metric to generate confinement and break confor- 
mal symmetry. In the meanwhile, the chiral symmetry 
is broken and a mass scale is introduced to simulate the 
effect. Based on AdS/CFT correspondence, the holo- 
graphic light-front QCD model yields a first order de- 
scription of some hadronic spectra. This model is quite 
appealing and promising, since it has established a pro- 
found relationship between super string theory and QCD 
in low energy limit. In this model, very few parameters 
(cutoff parameter Aqcd) are used to obtain the spectra 
for both mesons and baryons, such as n, p, and A, etc., 
which fit the experimental data well 10]. However, for 
the large body of mesons, only few of them are described 
properly and a large part of mesons are still left over. 
Besides, in its preset form the full spin interactions are 
not treated properly although it has potential to describe 
spin splittings. 

provides a convenient 



The light-front formalism If 



nonperturbative framework for the relativistic descrip- 
tion of hadrons in terms of quark and gluon de- 
grees of freedom[12l|. Some fundamental nonpertur- 
bative light-front QCD approaches are available, such 
as light- front Bethe-Salpeter approach [13j. holographic 
light-front QCD model[10], and light-front Hamiltonian 
method 14]. The light-front Bethe-Salpeter approach has 
been proposed by Kisslinger et al. to study pion form 
factor and the transition from non-perturbative to per- 
turbative QCD calculation of pion form factor. Like the 
B-S approach of instant form, the equation of motion for 
light-front B-S wave function should be solved together 
with Schwinger-Dyson equation for dressed quark prop- 
agator, vertex, and self-energy, and the model parame- 
ters include confining potential strengths, and others for 
parametrization of the BS Kernel and the running quark 
masses. An interesting conclusion drawn from the study 
of this approach is that the perturbative QCD calcula- 
tion works at the energy of 4-5 GeV, much lower than 
that explored previously by the instant form of QCD. 

The effective light-front QCD Hamiltonian theory pro- 
posed by Brodsky and Pauli [14| is an attempt to de- 
scribe the hadron structure as a bound constituent quark 
system in terms of Fock-space for the light-front wave- 
function. The effective Hamiltonian of the approach 
has been constructed recursively from the larger valence 
quark and anti-quark Fock sectors and reduced to the 
lowest valence quark- anti-quark sector [15j. Because of 
some unique features, particularly the apparent simplic- 



ity of the light-front vacuum, this model is a promis- 
ing approach to the bound-state problem of relativis- 
tic composite systems. Within the framework of the 
discretized light-front QCD, Pauli etal. have derived 
non-perturbatively an effective light-front Hamiltonian 
for mesons, which acts only on the qq sector [la, LL2J- 
The mass eigen equations of mesons are formulated in 
momentum-helicity representation which hinders its so- 
lution in total angular momentum representation. Be- 
sides, in this effective Hamiltonian, confining potentials 
and flavor mixing interactions are lacking, so that the 
excited states of mesons and flavor diagonal light mesons 
can not be treated properly [18[. 

In order to apply the approach to describe mesons in 
whole qq sector, essential changes are needed. First we 
have proved that in center of mass frame (rest frame) 
and in internal coordinate Hilbert subspace, the total 
angular momentum of the meson system is conserved(see 
Appendix A and B ). Then we are working in center of 
mass frame and in internal coordinate Hilbert subspace 
and make the following three significant improvements on 
the model: (1) transforming mass eigen equations from 
momentum-spin representation to total angular momen- 
tum representation and establishing a set of coupled ra- 
dial mass eigen equations for each total angular momen- 
tum ; (2) introducing a relativistic confining potential 
into the effective meson interaction phenomenologically 
based on lattice QCD results ; (3) including an SU(3) 
flavor-mixing interaction in the model phenomenologi- 
cally and obtaining a set of coupled radial eigen equa- 
tions for different flavor components. In having done 
above, finally we have a complete QCD inspired relativis- 
tic bound state model for mesons on the whole qq sector. 
This model has been applied to about 265 mesons with 
available data and with total angular momentum from 
J = to 6. The mass spectra, squared radii, and decay 
constants are calculated, and the calculated results are 
in good agreement with the data. While the most im- 
portant physical results have been reported briefly in a 
short letter [19], the present article will provide detailed 
information and solid foundation of the model for com- 
pletion 

This paper is organized as follows. In Sec. II the 
QCD inspired relativistic bound state model for mesons 
is described and the relativistic mass eigen equations for 
bound states with any total angular momentum are de- 
rived. In Sec. Ill based on lattice QCD results, a rela- 
tivistic confining potential in momentum space is intro- 
duced in the effective interaction of mesons. The effective 
interaction is extended to include an SU(3) flavor mixing 
interaction in Sec. IV. In Sec. V we present the numerical 
solutions for 265 mesons including both flavor-off and fla- 
vor diagonal mesons with J = — 6. Sec. VI is an analysis 
of the results obtained. Finally, conclusion and discus- 
sion are given in Sec. VII. The four Appendices are for 
clarifying some important issues and for the derivation 
of key equations. 



II. 



DESCRIPTION OF THE MODEL 



For convenience, Brodsky and Pauli defined a light- 
front Lorentz invariant Hamiltonian [14J 



H, 



LC 



P^P, 



p-p^ 



P? - Ml 



(1) 



The relativistic bound state problem in front form can 
be solved by solving the light-front mass eigen equation: 



# LC |*) 



M n 2 |*) 



(2) 



If one disregards possible zero modes and works in the 
light-front gauge, this equation can be solved in terms of 
a complete set of Fock states |// n ): 

]T f d[v! n ,]{n n \H hC \lj!nWn>\*) = M%{» n \*) . (3) 



For a meson, the ket |\&) holds: 

|*mcson) = ^ ^qq(xi, k±i, Xi)\qq~) 
i 

i 

+ ^2 ^qqg( x ^k±i, \i)\qqg) 



,^gqqq(xi,kxi,Xi)\qqqq) 



+ . 



(4) 



Within the framework of discrete quantization of light- 
front QCD, infinite dimensional Fock space has been 
truncated at a proper cutoff energy and the energy trun- 
cation plays a role of renormalization in discrete light- 
front QCD. By Tamm-Dancoff projection method and 
resolvent technique, the equation of motion in a larger 
Fock space of multi-particles can be reduced to that in 
a smaller one with an effective interaction to account for 
the effect of the projected out part of the Fock space. 
The reduction and projection procedure can be carried 
out recursively, finally the effective Hamiltonian and its 
eigen equation on qq sector can be obtained. For flavor 
off-diagonal mesons, disregarding the zero modes and the 
two-gluon annihilation effect, Pauli et al. has obtained 
the effective mass eigen equation for mesons in light-front 
relative momentum coordinate space [16J, LLZf 
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This is a relativistic covariant mass eigen equation for 
mesons in center of mass fame and in internal Hilbert 
subspace. However the equation of motion is written 
in relative momentum and helicity representation, and 
the momentum-helicity plane wave function contains all 
possible components of partial waves of the spin spherical 
harmonic functions $>ji s m , the total angular momentum 
J and its z-component M are not conserved. 

Despite this, Trittmann and Pauli[20] found an appro- 
priate method which can calculate the eigenvalue spec- 
trum separately for each J z = M. To do so, they trans- 
formed the light-front coordinate x back to the coor- 
dinate &3 by Terent'ev transformation [2 1| . and used a 



unitary transformation to transform the Lepage-Brodsky 
spinors to the Bjorkcn-Drell_spinors[22|- Then the mass 
eigen equation ([5]) becomes 23] : 



M 2 - (E t (k) + E 2 (k)f 
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This integration equation is written in momentum- 
spin representation in terms of internal relative momenta 
of two quarks, spin singlet and triplet are mixed. For 
the same reason as discussed above, the momentum-spin 
plane wave does not conserve J and M. As noted in 
Rcf[14], in general it is difficult to explicitly compute 
the total angular momentum of a bound state by us- 
ing light-front quantization. However, as addressed in 
Introduction, in the center of mass frame and in internal 
Hilbert subspace, the total angular momentum is con- 
served. This makes it possible to solve the mass eigen 
equation in total angular momentum representation (see 
Appendices B,C,D ). 

Since in center of mass frame and in internal Hilbert 
subspace, the total angular momentum J 2 and J z are 
conserved, we can transform the mass eigen equation 
(6) from momentum-spin representation to total angular 
momentum representation and establish the mass eigen 
equation for each J. Expanding the momentum-spin 
plane wave function in terms of the spin spherical har- 
monic functions ^jsiM^k, si, s 2 ) and projecting out the 
spin and angular part of the wave function in | JslM) 
subspace by the projecting operation, 

^2^2(lmsfi\JM)(-s 1 -s 2 \sn)Y lm (n k )x(s 1 ),x(s2) 



m/i S1S2 



(7) 

we obtain the mass eigen equation for the radial wave 
function of Rj s i(k) (see Appendix C). 



J+s' 



Rjal(k) 
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= E E J k' 2 dk'U J sl . slll {k;k')RjMk'). 



l' = \J-s'\ s'=0,l 



This is a set of coupled equations for radial functions 
Rjsi(k) of different partial waves and of spin singlet and 



triplets, coupled by the tensor potential and by the rel- 
ativistic spin-orbital potential. In this case, the eigen 
wave functions Rj s i(k) has the conventional definition 
and physical meaning. The bound states of mesons can 
be described concisely by the spectroscopic symbol of 
2S+1 Lj. 

The kernel Ug l . s , l ,(k;k') can be written as (see Ap- 
pendix D ), 

uL-v (*; k') = Y, E E / / dn " dn "' 

mm' S1S2 a^s^ 

x (Y lm (n k )\u aia2;sW2 (k,k')\Y llrn ,(n k ,)} (9) 

x (lmsfi\JM)(^s 1 ^s 2 \s f i)(l'm's'^\JM)(^s' 1 ^s' 2 \s'fi'). 

The above kernel U al . sn , (k; k') contains different kinds 
of central potentials, relativistic spin-orbit coupling po- 
tentials, and tensor potentials changing I by AZ = ±2 
and mixing spin singlet and triplets (see Appendix D ). 



III. INTRODUCING A CONFINING 
POTENTIAL 

Quark confinement is one of the fundamental problems 
in QCD for hadronic physics. The confinement and the 
spontaneous breaking of chiral symmetry are key ingre- 
dients for solving the low-energy hadronic bound states 
from QCD, but none of them has been completely un- 
derstood and solved. Numerical results show that the ef- 
fective light-front Hamiltonian model proposed by Pauli 
et al. without confining potentials can well describe the 
ground states but can not apply to the radial excited 
states of mesons. To describe the excited states properly, 
the confining potential must be included in the model|18|. 

Fortunately, we can refer to the constituent quark 
model which is successful due to the inclusion of a phe- 
nomenological confining potential in some way 24j. The 
key idea of this model consists in the introduction of a 
linear confining potential in coordinate space based on 
the numerical calculations of lattice QCD, and this non- 
relativistic confining potential can be generalized to rel- 
ativistic form. 

In nonrelativistic quark models the confining potential 
in configuration space is, 



V con {r) = \r + c 



(10) 



where A is the strength of the linear interaction, and c 
is a constant irrelevant in the present case and omitted 
hereafter. By Fourier transformation, the counterpart of 
the linear term Xr in momentum space is obtained, 



Min(g) - - 

q = k k! 



1 



(11) 



At the point of q — 0, the singularity indicates that the 
directly transformed result of linear potential could not 



be described correctly in momentum space, which results 
in an ill-defined bound state equation[25]. However, some 
different methods were employed to solve this problem 
for the relativistic case. In the present paper, the correct 
form for Vy m (q) is constructed by introducing a small 
parameter 77: 



V(q) = lim 

>)^o 2tt z or} 1 



1 



V 



(12) 



The relativistic linear potential in momentum space 
V\i n {Q) is a direct generalization of the nonrelativistic 
one, just replacing the nonrelativistic \q\ 2 in (|12[) by the 
relativistic Q 2 , which has the following specification in 
6. 261, 



Q 2 = (k - k' 



and 



w 2 = (E 1 -E[)(E 2 -E' 2 ) 
Then the form of relativistic confining potential is, 

1 
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(14) 



(15) 



Obviously, this confining potential is Lorentz covariant 

and can be used in either spin system or non-spin system. 

Now as the relativistic confining potential V con (Q) 

is included in the interaction, one has the new kernel 



U, 



,(k,k'), 



U Sl s 2 ;s' 1 s' 2 {' i ; k ) — 




hj\ Cj 2 ZVi £/ 2 



(16) 
-k, s 2 ) 



x [ 7 « . 7 ( 2 )"W + I* 1 ) • l^V s }u(k', s' 1 )u(-k', 4017) 
The scalar and vector interaction potentials read 
a(Q) 3 
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where e represents the scalar- vector mixing of the confin- 
ing potential. 



IV. INCLUDING A FLAVOR MIXING 
INTERACTION 

It is extremely difficult to derive a simple form of fla- 
vor mixing interaction in the above effective Hamiltonian 
from light-front QCD at present. However, without fla- 
vor mixing potential, one can not deal with the flavor 
diagonal mesons such as 7r°, p°, and / , etc. In the fun- 
damental hadronic theory, the quarks of u, d, and s have 
an approximate SU(3) symmetry. Due to this symmetry, 




the quarks fields transform each other under the SU(3) 
transformation [271 ]. 



exp[tJ2(0lT a + e^T al5 )} d 



where T a are Gell-Mann Matrices. For convenience of 
numerical calculation, we introduce phenomenologically 
a simple flavor mixing interaction as follows, 

Vf =7o [T+(l)T+(2) + T-(1)T"(2)] 
+<5 [T+(l)T+(2) +T-(l)T-(2) 
+ T+(l)T+(2) + T- s (l)T- s (2)] (19) 

where 70 and So are the strengths of flavor-mixing inter- 
action, the index 1 and 2 denote the quark and anti-quark 
in meson, respectively. The flavor SU(3) wave functions 
and generators are defined as, 



(20) 
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(22) 



(23) 



(24) 



The action of the flavor mixing interaction on flavor 
wave function is as follows, 



V f \uu) -- 


= 7o|dd) +(5o|ss) 


(25) 
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= §0 1 dd) + Sq | uu) 


(27) 



Combining this interaction with the precious one in 
equations (18) , we have a set of flavor-coupled radial 
eigen equations for the flavor components of up, down, 
and strange quarks, 



M 2 - (E^k) + E 2 (k)Y 



J+s' 



Rf s f{k) 



(28) 
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The interaction kernel including the flavor-mixing inter- 
action is 
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where W PlP2 '' p A 2 (K k') is defined as 






SiS2',S 1 S 2 




u(pi,k, *i)tt(p2, -k, s 2 ) [ 7 M (pi) • 7m(P2)W + ^ (1) ' ^ (2 V S 
[J, + V f ] u'(p[X, s'Ju'ip'z, -k', 4). (30) 



where piP2 1 p'iP , 2 — {uu,dd,ss}, If is the identity oper- 
ator in flavor space, Vy and Vs are the vector potential 
and scalar potential, respectively. For the flavor mixing 
mesons with total angular momentum J, the mass eigen 
equations are described explicitly by the following set of 
flavor-coupled equations 



Ml - (E u (k) + E a (k)) 2 \ Rf sl {k) (31a) 
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V. NUMERICAL SOLUTIONS 

In the above equations, J, s, Z denote total angular mo- 
mentum, total spin, and total orbital angular momen- 
tum, respectively. Different mesons can be classified by 
the spectroscopic symbol 2S+1 Lj(or their combination), 
which is equivalent to the symbol J PC . The space parity 
and charge conjugation parity are denoted as, 



P = (-1) L+1 
C = (-1) L+S . 



(32) 



In the present model, the mass eigen value problem of 
mesons is described by a set of coupled integration equa- 
tions, and the interaction includes a quark-anti-quark one 



gluon exchange potential Voce , a confining potential 
V con , and a flavor mixing interaction Vf. If the confin- 
ing potential has a pure iso-scalar structure( e = 0) and 
the flavor mixing interaction is omitted, the interaction 
contains two parameters: the effective coupling constant 
a and the confining potential strength A. Besides, the 
flavor mixing interaction has two parameters, and the 
constituent quark masses are also indispensable parame- 
ters to describe spontaneous chiral symmetry breaking. 
The numerical solution of the eigen equations can be 
x obtained by discretization of integration equation (8) or 
(31a — c), and the integration equations are transformed 
into matrix equations. The 4-fold integral of the kennel 
is completed by the integration technique of spherical 
harmonic functions and the angular momentum algebra, 
while the integration over k is performed by using Gauss- 
Legendre quadratures. The integration region k £ [0, oo) 
is projected onto the finite interval x G [—1,1] by x = 
j^j . The radial mass eigen equation (8) is discretized as 
follows: 



M 2 - (E 1 (k i ) + E 2 {ki)) 2 RjAk 



(33) 
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N 



= E E ^l U 'si;s>i>{ki\k j )Rjs>i>{k 1 )k 2 w, J , 

l' = \J-s'\ s'=0,l j=l 

where Wj is the weight of integration. Before diagonal- 
izing this matrix equations, special care should be taken 
for two kinds of singularities: the singularity at infinite k 
and the singularity as fc = k! inside the region of in- 
tegration. The first one has been solved by the pro- 
jection of the region k £ [0, oo) onto the finite interval 
x e [—1,1], and the second one is treated by infrared 
singularity treatment. The detailed procedures of calcu- 
lation can be found in Ref . [l8|, [28| . 

The parameters of the model are determined from best 
fit to experimental data. In this paper, a purely scalar 
confining potential(e = 0) is used. Reproducing the 
masses of 7r°, tt , and 7r(1300), we can determine a, A, 
and the masses of up and down quarks. Then by repro- 
ducing the masses of K , D°, and B , the mass pa- 
rameters of strange, charm, and bottom quarks are ob- 
tained. The parameters of flavor mixing interaction are 
determined by the best fit to the data of flavor diagonal 
mesons. From all the available data of mesons [23] with 
J = — 6 ( 12 mesons are left for future study: including 
6 exotic mesons and 6 mesons without any information 
about their J,s,L ), we have obtained an appropriate 
set of 6 parameters for flavor off-diagonal mesons: a = 
0.2574, A = 0.92 x 10 4 MeV 2 , m u/d = 0.297GeV,m s = 
0A18GeV,m c = lM3GeV,m b = i.U7GeV; and for the 
flavor diagonal mesons: 70 = 0.1 and So = 0.1. The 
number of the model parameters is minimum for this 
kind of semi-phenomenological models and comparable 
to BSE and CQM. The masses and wave functions of 
scalar and pseudoscalar, vector and axial- vector, tensor 
and pseudotensor mesons, and others with J = 3 — 6 
have been calculated and compared with the experimen- 



tal data in the Table( including 265 mesons and anti- 
mesons: 123 (u,d)-light mesons, 50 (s,u/d)-K mesons, 
24 (c,u/d)-D mesons, 14 (s,c)-D mesons, 12 (b,u/d)-B 
mesons, 10 (s,b)-B s mesons, 2 (c,b)-B c mesons, 16 (c,c) 
mesons, and 14 (b,6 ) mesons ). It is remarkable that 
among 265 mesons, 259 mesons are well described by 
this model within mass error less than 23%. 

In addition, the radial wave- functions of mesons in con- 
figuration space can be obtained from the radial wave- 
functions in momentum space by Fourier transforma- 
tion(see Appendix C), then one can calculate the mean 
square radii and the decay constants for some pseu- 
doscalar mesons listed in TabU and compared with ex- 
perimental data. 



TABLE I: The mean square radii and decay constants of 
some pseudoscalar mesons, compared with the experimental 
data[29|]. (Radii are given in fm 2 and decay constants are 
given in MeV) 





7T+ 


K+ 


D+ 


D s 


B 


(r 2 )the 


0.385 


0.253 


0.235 


- 


- 


y ) exp 


0.452 


0.314 


- 


- 


- 


fthe 


135.2 


210.7 


189.2 


253.1 


227.5 


J exp 


130.4 


155.5 


205.8 


273 


216 



TABLE II: The pseudoscalar mesons mass spectra (in MeV). 
~~Meson I U (J^ U ) Exp(Mev) Our , s(Mev) err(%) 



V 
?7(958) 
7/(1295) 

tt(1300)* 
tt(1300)° 
r/(1405) 
7/(1475) 
7/(1760) 

77(1800)* 

7r(1800)° 
X(1835) 

7/(2225) 

fo 
#(1460) 
#(1830) 

D° + 

D= 

Df 

B + 
7? C (2S) 



TABLE III: The scalar mesons mass spectra (in MeV). 
Meson J^° Exp(Mev) Our's(Mev) err(%) 



0"+ 
0"+ 


135 
140 

548 


135 
140 
143 





73 


o-+ 


958 


690 


27 


o-+ 


1294 


1258 


2.8 


o-+ 


1300±100 


1408 





o-+ 


1300±100 


1350 





o-+ 


1410 


1652 


17 


o-+ 


1476 


1700 


15 


o-+ 


1756 


1769 


0.8 


o-+ 


1816 


1454 


19 


o-+ 


1816 


2096 


15.4 


3-+? 


1833 


2110 


15.1 


o-+ 


2220 


2160 


2.7 


_ 
(T 
(T 


494 
498 
1460 


494 
494 
1522 




0.8 
4.2 


0" 


1830 


1597 


12.7 


<r 

0" 
0" 


1865 
1869 
1969 


1931 
1931 
2001 


3.5 
3.3 
1.6 


0" 

o - 


5279 
5279 
5367 


5584 
5584 
5667 


5.8 
5.8 
5.6 


<r 


6286 


6342 


0.9 


o-+ 


2980 


2980 





o-+ 


3637 


3533 


2.9 


o-+ 


9391 


8800 


6.3 



/o(600) 


0++ 


400-1200 


736 





/o(980) 


0++ 


980 


994 


1.4 


a (980)° 


0++ 


985 


1080 


9.6 


00(980)* 


0++ 


985 


930 


5.6 


/o(1370) 


0++ 


1200-1500 


1231 





o (1450)° 


0++ 


1474 


1333 


9.5 


00(1450)* 


0++ 


1474 


1457 


1.1 


/o(1500) 


0++ 


1505 


1522 


1.1 


/o(1710) 


0++ 


1724 


1568 


9.0 


/o(2020) 


0++ 


1992 


1606 


19.4 


/o(2100) 


0++ 


2103 


1989 


5.4 


/o(2200) 


0++ 


2189 


2026 


7.4 


/o(2330) 


0++ 


2321 


2052 


11.6 


# *(800) 


0+ 


672 


731 


8.8 


#o*(1430) 


0+ 


1412 


1535 


8.7 


Do (2400)° 


0+ 


2352 


2254 


4.2 


L>5(2400)* 


0+ 


2403 


2254 


6.2 


^0(2317)* 


0+ 


2317 


2169 


6.4 


Xco(LP) 


0++ 


3415 


3352 


1.8 


Xm(LP) 


0++ 


9860 


9860 





Xw(2P) 


0++ 


10232 


9990 


2.3 



TABLE IV: The axial vector meson mass spectra (in MeV). 



Meson 


J PU 


Exp(Mev) 


Our's(Mev) 


err(%) 


fa(1170) 


1+- 


1170 


1027 


12.2 


6i(1235)° 


1+- 


1230 


1127 


8.4 


oi(1235)* 


1+- 


1229 


1343 


9.3 


ai(1260)° 


1++ 


1230 


1276 


3.7 


ai(1260)* 


1++ 


1230 


1371 


11.4 


A (1285) 


1++ 


1281 


1295 


1.1 


fa (1380) 


1+- 


1386 


1301 


6.1 


/1 (1420) 


1++ 


1426 


1311 


8.0 


/1 (1510) 


1++ 


1518 


1419 


6.5 


fa (1595) 


1+- 


1594 


1495 


6.2 


01 (1640)° 


1++ 


1647 


1745 


6.0 


ai(1640)* 


1++ 


1647 


1724 


4.7 


#i(1270) 


1+ 


1273 


1459 


14.6 


#i(1400) 


1+ 


1402 


1484 


5.8 


#i(1650) 


1+ 


1650 


1757 


6.4 


Di(2420)° 


1+ 


2422 


2400 


0.9 


Di(2420)* 


1+? 


2423 


2400 


0.9 


£>i(2430)° 


1+ 


2427 


2425 


0.1 


Dsi(2460)* 


1+ 


2460 


2530 


2.9 


L>si(2536)* 


1+ 


2535 


2549 


0.6 


Bi(5721)° 


1+ 


5721 


5666 


1.0 


Bsi(5830)° 


1+ 


5829 


5800 


0.5 


Xoi(lp) 


1++ 


3510 


3504 


0.2 


fai(lp) 


1+- 


3526 


3509 


0.5 


XbiO-p) 


1++ 


9892 


10040 


1.5 


Xbi(2p) 


1++ 


10255 


10040 


2.1 



TABLE V: The vector meson mass spectra (in MeV). 



Meson 


J PC 


Exp(Mev) 


Our's(Mev) 


err(%) 


p(770) u 


1~ 


775 


1015 


31 


P (770)± 


1— 


775 


1239 


60 


w(782) 


1" 


783 


1270 


62 


<^(1020) 


1" 


1019 


1334 


31 


oj(1420) 


l" 


1425 


1410 


1.0 


(0(1450)° 


1— 


1465 


1636 


11.6 


,9(1450)* 


1~~ 


1465 


1323 


9.7 


p(1570)° 


I— 


1570 


1641 


4.5 


0(1570)* 


1~~ 


1570 


1740 


10.8 



w(1650) 


1"" 


1670 


1675 


0.3 


0(1680) 


1— 


1680 


1786 


6.3 


p(1700)° 


1— 


1720 


1836 


6.7 


p(1700)* 


1"" 


1700 


1362 


19.8 


p(1900)° 


1— 


1909 


1996 


4.6 


p(1900)* 


1" 


1909 


1761 


7.8 


,0(2150)° 


1~ 


2149 


2087 


2.9 


p(2150)* 


1 __ 


2149 


2430 


13.1 


K* (892) 


1~ 


892 


1345 


50.1 


K*(UW) 


1" 


1414 


1415 


0.1 


#*(1630) 


1~? 


1629 


1502 


7.8 


#* (1680) 


1" 


1717 


1531 


10.8 


73* (2007)° 


1" 


2007 


2100 


4.6 


73*(2010)* 


1" 


2010 


2126 


5.8 


73* (2640) 


1"? 


2637 


2403 


8.9 


Df 


1"? 


2112 


2214 


4.8 


D S i(2700) ± 


1" 


2690 


2233 


16.9 


B* 


1" 


5325 


5518 


3.6 


£* 


1" 


5413 


5625 


3.9 


J/V 


1— 


3097 


3284 


6.0 


V(2S) 


1"" 


3686 


3362 


8.8 


■0(3770) 


1~ 


3773 


3684 


2.3 


V>(4040) 


1— 


4039 


3700 


8.4 


V>(4160) 


1— 


4153 


4197 


1.1 


X(4260) 


1— 


4263 


4769 


11.9 


#(4360) 


l~~ 


4361 


4783 


9.7 


^(4415) 


I— 


4421 


5341 


20.8 


#(4660) 


l~~ 


4664 


5418 


16.2 


7(15) 


1— 


9460 


9693 


2.4 


7(25) 


1~~ 


10023 


10013 


0.1 


7(35) 


I— 


10355 


10060 


2.8 


7(45) 


I— 


10580 


10765 


1.7 


7(10860) 


1~~ 


10865 


10783 


0.8 


7(11020) 


l~~ 


11019 


10861 


1.4 



TABLE VI: The tensor and pseudotensor meson mass (in MeV). 



Meson 


J PU 


Exp(MeV) Theor(MeV) 


err(%) 


tt 2 (1670) 


2"+ 


1672 


1587 


5.1 


tt 2 (1880) 


2" + 


1895 


1589 


16.1 


7T 2 (2100) 


2" + 


2090 


1922 


8.0 


#2 (1580) 


2~ 


1580 


1530 


3.2 


#2(1770) 


2 - 


1773 


1539 


13.2 


#2(1820) 


2" 


1816 


1763 


2.9 


# 2 (2250) 


2 - 


2247 


1765 


21.4 


7(173) 


2" 


10161 


10218 


0.6 


a 2 (1320) 


2++ 


1318 


1421 


7.8 


a 2 (1700) 


2++ 


1723 


1474 


14.5 


#2* (1430) 


2+ 


1425 


1531 


7.4 


#2(1980) 


2+ 


1973 


1575 


20.1 


731(2460)* 


2+ 


2460 


2456 


0.1 


732(2460)° 


2+ 


2462 


2456 


0.2 


73^(2573)* 


2+? 


2573 


2580 


0.3 


73} (5732) 


2+? 


5698 


5706 


0.1 


55(5747)° 


2+ 


5743 


5765 


0.4 


BJ 2 (5840)° 


2+ 


5840 


5831 


0.2 


5^(5850)° 


2+? 


5853 


5883 


0.5 


Xc2(lP) 


2++ 


3556 


3732 


4.9 


Xc2(2P) 


2++ 


3929 


3745 


4.7 


Xa(lP) 


2++ 


9912 


10014 


1.0 


Xb2(2P) 


2++ 


10269 


10354 


0.8 



^ 3 (1670) 
p 3 (1690)* 
p 3 (1690)° 
03(1850) 
p 3 (1990)* 
p 3 (1990)° 
p 3 (2250)* 
p 3 (2250)° 
# 3 *(1780) 
# 3 (2320) 
04(2040)* 
a 4 (2040)° 
/ 4 (2050) 
/ 4 (2300) 
#1(2045) 
# 4 (2500) 
,95(2350)* 
ps(2350)° 
# 5 * (2380) 
a 6 (2450)* 
a 6 (2450)° 
/ 6 (2510) 



3— 


1672 


1677 


0.3 


3— 


1688 


1702 


0.8 


3— 


1688 


1697 


0.5 


3— 


1854 


1807 


2.5 


3— 


1982 


1795 


9.4 


3— 


1982 


1807 


8.8 


3— 


2230 


2660 


19.2 


3~~ 


2230 


1852 


17.0 


3- 


2324 


1777 


23.5 


3+ 


2324 


1812 


22.0 


4++ 


2001 


1745 


12.7 


4++ 


2001 


1743 


12.9 


4++ 


2018 


1865 


7.6 


4++ 


2300 


2016 


12.3 


4+ 


2045 


1827 


10.6 


4" 


2490 


1933 


22.3 


5~~ 


2330 


2292 


1.6 


5"" 


2330 


2218 


4.8 


5" 


2382 


2352 


1.3 


6++ 


2450 ± 130 


2412 





6++ 


2450 ± 130 


2423 





6++ 


2465 ± 50 


2649 


5.3 



TABLE VII: The mesons of J > 3 (in MeV). 
Meson J ru Exp(MeV) Theor(MeV) err(%) 



VI. ANALYSIS OF THE RESULTS 

In the above calculations, only one set of parameters 
are used, which deserves discussion. The effective cou- 
pling strength or running coupling constant a and the 
related constituent quark masses have a great influence 
on ground state of light mesons, such as n. The confin- 
ing potential strength A governs the quark confinement 
at large distances and has strong influence on the ex- 
cited states of light mesons and also on the spectra of 
heavy mesons. From the recent experiments of hadron 
physics, we know that the QCD coupling a(Q 2 ) becomes 
large constant(not singular) in the low momentum limit, 
which is called infrared conformal invariance 30]. This 
experimental fact explains why our model with a set of 
constant parameters works well to describe the structures 
of mesons in the energy region of 0.14GeV— > lOGcV, and 
our results may be thought of confirming the infrared 
conformal invariance feature of QCD on meson sector. 

For light scalar mesons such as ao ,Kq, etc., although 
the structure of the scalar mesons remains a challeng- 
ing puzzle, our model still describes ao(980), a (1450), 
#q(800), etc. quite well. For heavy mesons, be- 
cause of the large masses of heavy quarks, the effec- 
tive double-gluon-exchange interactions for off-diagonal 
heavy mesons are weak, which makes the model applica- 
ble to them. Therefore, the calculated mass spectra for 
the mesons of u/ds, u/dc, sc, cc, cb, u/db, sb, and bb are 
in good agreement with the data. However the meson 
#*(892) on u/ds sector with larger error of 50.1% needs 
special investigation(see below). 

It should be noted that the J and P of D* are not 
identified by experiments, but their width and decay 
modes are observed and consistent with the 1~ state. 
Nevertheless, our model provides a definite assignment 



of J = 1 and P = — 1 for D*^. A similar prediction of 
the unidentified J and P is also made for other 8 mesons: 
A(1835), Z?i(2420) ± , IT* (1630), L>*(2640), J DJ 2 (2573) ± , 
B}(5732), 5^(5850)°, and /j(2220). 

The 6 mesons with errors larger than 23% provide some 
information. For the vector mesons of rj, 77' (985), p(770)°, 
0(1020), and w(782) on u/d sector, and if* (892) on 
(u/d)s sector , the large discrepancy indicates that the 
structures of these mesons are special than others and 
need a different set of parameters: indeed, as the set of 
parameters are re-adjusted to the set of (a — 0.4594, 7 = 
0.58, S — 0.74) and with the others the same, a better 
fit is found with errors less than 23%. Increase of the 
effective interaction strengths implies that these vector 
mesons may have strong coupling between qq and qqqq 
subspaces and among different flavor components. 



VII. CONCLUSION AND DISCUSSION 



agreement with the data. 



2. Comparing to other BSE and CQM meson models, 
our model is novel in following points: 1) The effective 
Hamiltonian is derived within the framework of light- 
front QCD and the form (the spinor structure ) of the 
effective interactions is fixed by the lowest order of light- 
front QCD. 2) The mass eigen equation is for the squared 
rest mass, the separation between kinematical energy op- 
erator and interaction operators is rigorous. 3)The spinor 
structure of the effective interaction make it momentum- 
energy dependent. 4) Also due to the spinor structure 
of the effective interactions, the dynamics of spin-spin, 
spin-orbital, and tensor interactions ( especially the spin 
singlet-triplet mxing and orbital angular momentum mix- 
ing ) are included ( see Appendix C,D ). 5) The predic- 
tive power and the descriptive precision of the model are 
much better. 



In conclusion, we have formulated the QCD inspired 
relativistic bound state model for mesons and derived its 
mass eigen equations in total angular momentum repre- 
sentation. It is proved that in center of mass frame and 
in internal Hilbert subspace, total angular momentum 
of the meson system is conserved. Moreover, by taking 
the advantages of other effective QCD approaches |a.[26|. 
the model has been improved significantly by introducing 
both a relativistic confining potential and an SU(3) flavor 
mixing interaction. The resulting radial mass eigen equa- 
tions are solved numerically and nonperturbatively, and 
265 mesons including flavor off-diagonal mesons and fla- 
vor diagonal ones with J = — 6 are calculated and com- 
pared with experimental data. The calculated masses are 
in good agreement with the data within the mean square 
root mass error of 14%, only 6 mesons with mass error 
larger than 23%. Besides, the wave functions obtained 
from the model also yield reasonable mean square radii 
and decay constants for some pseudo scalar mesons. In 
view that the structure of the light scalar mesons is still a 
subject of controversy [3lj, and the internal dynamics of 
heavy-light mesons in the static limit is far more compli- 
cated than that of the heavy-heavy ones [3^. our model 
can be thought to be successful to describe a large body 
of mesons. 

The comparison of our model with other approaches is 
as follows: 

1. As Pauli's model is concerned, we have improved the 
model significantly on 5 important points and make it a 
predictive and systematic model for mesons: 1) Proving 
that in internal Hilbert subspace, total angular momen- 
tum is conserved; 2) establishing the mass eigen equa- 
tions in total angular representation for the first time; 3) 
introducing the relativistic confining potential into the 
model, which is new and quite different from Pauli, and 
its form taken from the [y, l26j |: 4) including the flavor 
mixing interaction; 5) solving the mass equations for 265 
mesons nonperturbatively and the results are in good 



3. Comparing to holographic light-font QCD model 
of Brodsky et q?.[10l|. our model has the following new 
aspects : 1) In the effective Hamiltonian of mesons, the 
kinematical energy operator is identical for both holo- 
graphic light-front QCD model and our model, but the in- 
teraction terms are quite different. 2) Holographic light- 
front QCD model does not specify the effective interac- 
tion in detail, but just simulates confining potential by 
boundary condition ( or harmonic oscillator potential), 
or recently by a positive-sign dilaton metric to generate 
confinement and break conformal symmetry; instead, our 
model provides a detailed semi-phenomenological effec- 
tive interaction including its spinor structure, the confin- 
ing potential, and the flavor mixing interactions. 3) Holo- 
graphic light-front QCD model does not include spin- 
spin, spin-orbital, and tensor interactions, the total an- 
gular momentum of the system is not treated properly ( 
although it has potential to describe the spin splittings 
); in the contrary, our model specifies the spin interac- 
tions and the spin dynamics is described fully in total an- 
gular momentum representation. 4) Finally, our model 
has been applied to a larger number of mesons ( 265 
mesons identified experimentally ) with higher precision 
than those of holographic light-front QCD model. In the 
above respects, our model has provided a tentative and 
effective solution to the problems listed above and the re- 
sults are amazingly in good agreement with experimental 
data. In this sense, our model can be considered to be 
of complementarity to and refinement of the holographic 
light-front QCD model. 
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Appendix A: Dynamics in light front form and 

instant form in center of mass frame and in internal 

Hilbert subspace 

To avoid misunderstanding of light-front dynamics, we 
start from a discussion of full contents of dynamics for 
both instant form(IF) and light front form(LF). The con- 
tent of dynamics should contain the following four as- 
pects, we list them for both dynamics of instant form 
and dynamics of light form as follows. 



1. Full contents of dynamics in instant form 

1) Definition of time x° = ct = t(c = 1) 

2) Hamiltonian (energy) operator is defined as the time 
translation operator: 



d 

ih Trr: ~ P° = H = M 
ox v 



(Al) 



M is dynamical mass operator. 

3) Dynamics 

(i) Time evolution dynamics: equation of motion 
(Schrodingcr equation), 



ih—- = HV = M* 
ox° 



(A2) 



(ii) Stationary dynamics: for stationary solution , 

tf (i) = e - iMt/h $ (A3) 

one has the Hamiltonian eigen equation 

H^> = Af* = M*, (A4) 

where M is the eigen value of M . 

4) Specification of dynamical operators and kinemati- 
cal operators among Poincare generators: 6 kincmatical 
operators: P l , J l ,(i = 1,2,3); 4 dynamical operators: 
P°,K\ (i = 1,2, 3). 

It should be noted that the dynamical operators con- 
tain interactions via the Hamiltonian and Lorentz boost 
operators while the kinematical operators do not. Conse- 
quently, the kinematical operators can be used to charac- 
terize the state of the system as good quantum numbers 
according their algebraic structure and the dynamical 
operators except the Hamiltonian operator can not play 
such a role. It should be emphasized that the above spec- 
ification is made in whole Hilbert space of the states of 
composite systems. For a composite many-body system, 
the whole Hilbert space of states can be factorized into 
two parts: a) the center of mass motion characterized by 
its momentum P , and b) the internal motion character- 
ized by internal quantum numbers and (J, J 3 ). Corre- 
spondingly, the Poincare operators contain two kinds of 
operations, one on the subspace of center of mass motion 



and the other on the subspace of internal motion. Since 
the center of mass motion can always be separated from 
the internal motion, the state wave function of the com- 
posite system ^ can be written as ^ = ^cm^ inter, where 
the wave function of center of mass motion is character- 
ized by the center of mass momentum, namely ^ C m = ^ p 
with P l \& p — P l ^> p , while the internal wave function is 
characterized by internal quantum numbers and (</, J 3 ). 



2. Full contents of dynamics in light front form: 

1) Definition of time x + : x + = x° + x 3 

2) Hamiltonian ("energy") operator is defined as the 
time translation operator: 



ih- 



d 

dx + 



(A5) 



From 



P + P- -Pl = PqP° - Pi - A 2 - H = Ma (A6) 



Mq is rest mass operator; one has 



(A7) 



3) Dynamics 

(i) Time evolution dynamics: equation of motion 
(Schrodingcr equation), 



*^ = -L ( m o2 + p 1) * 



(A8) 



(ii) Stationary dynamics: for stationary solution 
value of P~ , P + and P± are eigen values of P + , P±), 



^(j^t,P + ,Pj_,x + ) with quantum numbers: "energy" 
E~ = |qr and momentum P = (P + , Pj_)(E~ is the eigen 



*(*+) = e- M2 * + / p+h *(^ T ,P+,P 1 



(A9) 



One has mass eigen equation 

1 



,M 2 



M 2 M 2 -. , 

— *(— ,P+,P ± ) = — (M 2 +Pl)v&(^ T ,P+,^)10) 



M 2 -> - M 2 

M^(p^,P + , Pl) = (M 2 - p 2 )*(_ p+, PJA11) 

4) Specification of dynamical operators and kinemati- 
cal operators among Poincare generators: 7 kinematical 
operators: P+, J 3 , P l {i = 1,2),K 3 , E 1 = K 1 + J 2 , E 2 = 
K 2 — J 1 ; 3 dynamical operators: P~ , F 1 = K 1 — J 2 , 
F 2 = K 2 + J i . 
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3. Dynamics in center of mass frame and in 
internal Hilbert subspace for both forms of dynamics 



where To and Tq are constant shifts of proper time. One 
can choose the start point of time such that 



The internal structure of a composite system should 
be described in the rest frame as well as in the cor- 
responding internal Hilbert subspace. Since the center 
of mass frame always follows the center of mass motion 
of the system and the position of the center of mass of 
the system is at the origin of the frame, the wave func- 
tion of center of mass motion of the system should be 
* p_ , and the center of mass momentum and the center 
of mass coordinates of the system should be zero, namely 

< *P =0 l^i*p=o >= ° and < *p=ol f i*p=o >= °- In 
the center of mass frame, the Hilbert subspace of cen- 
ter of mass motion is frozen to ^ p =0 , the whole Hilbert 
space of states of the system is thus projected onto the 
corresponding internal Hilbert subspace Winter- Conse- 
quently, the dynamics of the composite system is reduced 
to the internal dynamics. Projecting onto the frozen cen- 
ter of mass wave function and integrating out the center 
of mass degrees of freedom, one obtain the Poincare op- 
erators in the internal subspace Winter as follows. 

1) Four momentum and property of time in center of 
mass frame and in internal Hilbert subspace. 

In center of mass frame, the wave function of center of 
mass motion is : ^ p =0 - The four momentum operator 
in internal Hilbert subspace can be obtained by project- 
ing out the center of mass degrees of freedom ( namely 
averaging over the center of mass wave function). Since 

P° = H = (M 2 + P 2 )V2 an d P l ^ p=Q = 0, in internal 
Hilbert subspace, one has 



p* = ( v&- \P l hjr- 

± inter \ x P=0 1 I P=0 



= 



Pt*„ = < *3 _JP |*5_n > = M 



P=0I ± l*P=0 



(A12) 



(A13) 



Here Mo is the operator of rest mass of the system. 
Thus in internal Hilbert subspace, the four momentum 
operators for instant form read: 

P? nter = (Mo, 0,0,0) (A14) 

while four momentum operators for light front form are: 

(M , 0,0, M ), (A15) 

From the above results , one has 

^inter = "inter = ^inter = M-Q, (AiOj 



inter 



This leads to 



and 



r = ->• x + = x v = 



TO = T {) = 



X = X = T 



(A19) 



(A20) 



(A21) 



2) Dynamics in center of mass frame and in internal 
Hilbert subspace 

(i) Time evolution dynamics: equations of motion in 
center of mass frame and in internal Hilbert subspace 

The Schrodingcr equations 



dx+ 

in light front form, and 
5* 



5* 

ift^r^- = M * 



m dx^* = Mo * 



in instant form become the same 



ih— = Mo* 

OT 



(A22) 



(A23) 



(A24) 



(ii) Stationary dynamics: mass (energy) eigen equa- 
tions in center of mass frame and in internal Hilbert sub- 
space. 

The mass eigen equations 

M * = Af * (A25) 

in instant form where Mo is the eigen value of Mo, and 

(A26) 



M 2 * = M 2 * 



in light front form are also the same because multiplying 
M on the first equation leads to the second one. 

3) Kinematical and dynamical operators in center of 
mass frame and in internal Hilbert subspace 

Projecting onto internal Hilbert subspace, the kine- 
matical and dynamical operators can be obtained from 
the following calculation. From the results of (A12-A15) 
and 



,J "cm ' "inters 



J\ m = x 2 P a - P 2 x 3 ,cycli 



(A27) 
(A28) 



ih— = ih— = h— 
dx + dx a dr ' 



(A17) 



where t is the proper time corresponding to the rest mass 
operator Mq. The last equation leads to 



X + = T + T , X° 



(A18) 



one obtain 



*p=ol <™l f=o/ — i). 



*P=ol J4 |*P=o) = J lnter 



(A29) 



(A30) 
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*P=ol^ 



*P=0 



*P=ol 



j,0p 4 _ $ tp 



(*P =0 l*i*P=o>^o = 0, 



w p= 



ol^l*p% 



w p= 



IX 1 + J 2 



*p= 



°l*pJ^3il) 



JirU&A^) 



manner and the results are the same as that in Appendix 
A. 

A relativistic dynamical system has inhomogeneous 
Lorcntz symmetry defined by the Poincare algebra: P^ 
is energy-momentum vector, and M^ v is used to de- 
scribes the rotational and boost transformations. In in- 
stant form, the angular momentum and boost vectors are 
given as: M* = e ijk J k and M oi = K l 

Now define the "quasi angular momentum" operators 
in the light-front form: 



*P=ol^ 



*P=0 



= <*P=o\K' 



-J 1 \y p=a ) = -J! 1 Wa) 



vp _ IF 1 !*- 



= (*P = ol^ 1 -^l 



^P=0 



* p= 



ol^l*p= 



*p= 



I if 2 



J 1 ^ !^^ 



From the above results, one obtains the same reduced 
and degenerated kinematical and dynamical operators for 
both forms of dynamics in internal Hilbert subspace as 
follows: kinematical operators: Ji n t er (i — 1,2,3); dy- 
namical operator: Pf nter — Pinter — Mo- 

The above results tell that in center of mass frame and 
in internal Hilbert subspace, light front time and instant 
time, light front dynamics and instant dynamics, light 
front angular momentum and instant angular momentum 
are identical. 



4. Conclusion 

In general frames and in whole Hilbert space, both 
forms of dynamics are quite different. However, in center 
of mass frame and in internal Hilbert subspace, the two 
forms of dynamics are reduced to the identical internal 
dynamics. 

There is a dilemma in this paper at first glance: our 
model begins with a light front QCD model, but the fi- 
nal form of our model possesses the feature of instant 
dynamics of QCD. Is it of LF dynamics or IF dynamics? 
The solution to the dilemma is given in this Appendix, 
the answer is that in center of mass frame and in internal 
Hilbert subspace, the reduced internal dynamics of both 
forms are identical. 

Therefore, our model contains ingredients of both the 
instant form and light front form of QCD, it can be called 
as QCD inspired effective Hamiltonian meson model. 



J 3 = J 3 



P+ 



'3r« 



-JSJ&34) 



Jirb@l35) 



J ±l = M - L e ij (-(F 3 ± P+ - E J ± P~) - K A P- 

+J 3 s 3l P[),(tJ = 1,2). 



(Bl) 



It is easy to prove that they satisfy the SU(2) algebra: 
[J i ,J j ]=ie ijk J k (B2) 

It is very useful to define a 'light-front Hamiltonian' as 
the operator: 



H LC = P»P U = P-P+-P 1 _ = M 2 



(B3) 



Hlc commutes with the quasi angular momentum op- 
erators : 



[H LC ,J}=Q. 



(B4) 



In principle, one could label the eigen states as 
\M, P + ,P±, J 2 ,Jz), since J3 is kinematical. However, 
J71 is dynamical and depends on the interactions. Thus 
it is generally difficult to explicitly compute the total spin 
J of a state using light-front quantization. Fortunately, 
in center-of-mass frame and in internal Hilbert subspace, 
by using the results of Appendix A, one has the following 
equations, 



(* P=0 \J 3 \*p= ) = J, 
(*p =0 \J i \*p =0 ) = J l 



3 

inter ? 



(B5) 



(M = l,2) 



Therefore, in internal Hilbert subspace, the quasi an- 
gular momentum operators Jl nter are identical to the 
total angular momentum operators Ji nter (i = 1,2,3), 
the total angular momentum is conserved, and the eigen 
equation of the Hamiltonian Hlc of the internal dynam- 
ics can be solved in the total angular momentum repre- 
sentation. 



Appendix B: Conservation of total angular 
momentum in internal Hilbert subspace 

The reduction of angular momentum operators in in- 
ternal Hilbert subspace can be discussed in an alternative 



Appendix C: Derivation of the radial mass eigen 
equations in total angular momentum representation 

According to Pauli et al., the effective mass eigen equa- 
tion of mesons of light-front QCD in center of mass frame 
and in internal Hilbert subspace reads: 
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M 2 - (E^k) + E 2 (k)) 2 VsiS2 {k) 



]T J d 3 kU SlS2;s{s , 2 (ft; k')v s[ s> 2 (*'). (CI) 



where 



U. 



4m s a(Q) 



3tt 2 Q 



R(Q) 



^siS2;s' 1 s' 2 
y/A(k)A(k') 



(02) 



with 



and 



S'aiaaiaisi = [ M ( fc > s l)7p( 1 )w(fc', s'l)] 

x {v(-k,s 2 )^(2)v(-k',s' 2 )} (C3) 



1 



1 



1 



A{k) ' r \Ek{k) E 2 (k)J' 

m\vn 2 



m s = 77i i + m 2 , m r 
Q = Q(k;k'). 



mi + m 2 ' 



(CM) 



Equation (Al) can be written as Schrbdinger equation 
in the light front QCD, 



77 ^ meson -^"0 ^nicson 



(C5) 



The general eigen wave function "JVeson of meson can be 
expressed in momentum-spin representation, 



*, 



Y J d 3 k VsiS2 (k) \ X (si)x(s2) ■ k) . (C6) 



Here basis of the momentum-spin representation are 



(r\x(si)x(s2) -k) = 



1 



(27T/7) 3 / 2 



X(si)x(s 2 )e 4fc - r , (C7) 



where the spin wave functions and their orthogonal con- 
ditions read 



<x(«i)lx(«» -«.,«• (C9) 

The orthogonal conditions of the spinors are 

(k ■ u(k, si)v(—k, s 2 )\u(k' , s[)v(—k', s' 2 ) ■ k') 

= S^(k-k')S Sls/i 5 S2S , 2 , (CIO) 

u(k,si)u(k,s[) = <5 SlS ;, (Gil) 

v(-k,s 2 )v(-k,s' 2 ) = 5 S2S , 2 , (C12) 

and the completeness conditions read, 

J2u(k,s)u(k,s) = J-( 7/ ,fcj*+ m ), (C13) 



J2v(-k,s)v(-k,s) = J-( 7/J fc£- m ), (C14) 



where fc£ = (£i(fc),fc), and fc£ = (E 2 (k),-k). 

According to the Dirac form of quantum mechanics, in 
the eigen equation (|C5[) . the Dirac form of the Hamilto- 
nian operator is 



H = E + U, 



(C15) 



where 



E = / d?k[E x {k) + E 2 {k)} 2 

J2\x(si)xM-k)(k- X (si)x(s2%(C16) 



and 



U= f d 3 kd 3 k' Y U(k,k') 



x [u(k, siM-fe, S2 )(^(i)rmHk', s[) V (-k', s y] 

x|x(si)x(s2)-fe)<fe'-x(si)x(4)l, (C17) 

with the definition, 

1 



*<>•*>- d?iP m 7mm 



(CIS) 



In the above equation, as done by Pauli et al.[20j. 
the light front k— space has been transformed back to 
the Lab k— space by the Terent'ev transformation, and 
Lepage-Brodsky (helicity) spinors have been transformed 
to the Bjorken-Drell (spin) spinors. 

Using eqs. (|C6l[C15IIC18| ) and projecting equation (|C5[) 
onto the subspace |x( s i)x( s 2) ■ k), we recover the equa- 
tion (jCip . indicating that the Dirac Form of the eigen 
equation (|C5|) is equivalent that of (|C1|) . 

Since (E±(k) + E 2 (k)) 2 and the interaction kernal op- 
erator U[k, k'\ <x(l),<r(2)] are scalar (see Appendix D, 
discussion below eq. (|D9jl ). H is rotational invariant with 
respect to the total angular momentum Ji = li+sj+sf = 
h + Si, [H,Ji] — 0. That means the total angular mo- 
mentum J and J z are conserved. Based on this point, 
the wave function of the meson system can be written in 
total angular representation as follows, 

_ ,7+s _ 
^meson(k, fife, s) = ^ Y ^ Rjsl{k)$ JslM (fife, s), 
J.M l=\J-s\ s=0,l 

(C19) 
were the total angular momentum eigen functions <&j s im 
of {J 2 , J z ,s 2 ,l 2 } are, 

$j s iM(ttk,s) = ^(?ms^|JM)F im (fife)xa M (12), (C20) 

mfj, 

the eigen wave functions of spin singlet and triplet read 

as. 

Xs/1 (12) = Y^ s ^ S2 \ s ^(si)x(s2). (C21) 
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By virtue of the Fourier transformation in spherical 
coordinates, from the eigen wave function in the momen- 
tum radial k— space, one can obtain the corresponding 
wave function in the configuration radial r— space, 

^M(r,fi„s) = Jdk 3 * JM (k,n k , s y kr 

■J+s „ 

= J2 E E / k 2 dkRj sl {k).h{kr) 
l=\J-s\ s=0,l V ,m' J 

dnfc*j.iiif(n fc ,*)y,f m ,(n fc )^ m /(n r ) 

J+s 

X] E iW)^Af(0-,s),(C22) 

i=|J-s| 8=0,1 



where 



$JsiM(Or,s) = ^(lmSLl\JM)Yi m (£lr)Xs^)i 



in ft 



R.Jsi(r) 



k 2 dkR., s i{k)Ji{kr) 



J(kr) = 1/442/ + 1) i l ji(kr). 



(C23) 



ji(kr) is the spherical Bessel function of order /. 

Using the expression (|C19| of the wave function 
^ meson, projecting the mass eigen equation (|C5j) onto 
the QjsiM subspace from the left, and integrating out the 
spin and angular part of the wave function, we obtain the 
eigen equations for the radial wave functions Rj s i{k), 



M£ - (E^k) + E 2 {k)Y R Jsl (k) 



J+s' 



(C24) 

E E / 'k a dk'u J al , alll {k;k')RjMk'), 

!'=|J-s'| a'=0,l J 

where the kernel Ugi. sil ,{k\ k') is defined as, 

U J sl . a , v (fc; fc') = E E E / / d ^ dn "' 

mm' s±s 2 a' x a' 2 

x (Yi m {n h )\U aiaa . a > 1 ^(k,Ts!)\Y Vm >(n k >)) (C25) 

x (lmsti\JM)(ls 1 ^s 2 \s^)(l'm's'fi'\JM)(^s , 1 ^s' 2 \s'^). 

This is a set of coupled equations for the radial functions 
Rjsi(k) that have different partial waves, spin singlet and 
triplet coupled by the tensor potentials and by the rela- 
tivistic spin-orbital potential (see below). 



u(k,s = — i) = 



v(—k, s 



v(—k, s 



where 



y/2m 


l(J01+TOi) 




1 


y/2m 2 (E 2 +m 2 ) 




1 



y/2m 2 (E 2 +m 2 ) 



( \ 

E\ + m\ 

V -A. I 

( ~ k * \ 
-h 

E 2 + m 2 

V 1 

( -fcr \ 

k z 



\ E 2 + m 2 / 



ki yr = k x ±ik y = ksm6 k e ±t ' Pk = kJ —Y 1±1 (9 k ,(pk), 

(Dl) 



/47T 



k z — kcoadk = k\j — Y w (0k, Vfe) 



Defining the spherical spinors 



*i.(fifc) = 53(lm^|| S >yoo(n fc )xM = -t=x(«), 



47T 



(D2) 



af.(n*) = V(imHi s )r lm (ft fc ) x O,) = -4=o*x(«), 

2 V 47T 



(D3) 
where a k = [cr ■ k)/k and f2j- = {Ok,(pk){<Jk is pseudo 
scalar ), the spinors can be re-expressed as 



u(k, s) 
v(—k,s) 



i4i(fc)*J B (n fc ) 

Bitfc)*!^) 

-B 2 (*)*| a (n fc ) 

A 2 (fc)$fffi fc ) 



(D4) 

(D5) 



Appendix D: Calculation of the interaction kernel in 
total angular momentum representation 

The quark and anti-quark spinors are given in the 
Bj0rken-Drell representation, 

( El + TOl \ 



k z 

V h ) 



u(k, s 



y /2m 1 (E 1 +mi) 



where 



Mk) = j MgLtl*) , Bi(fc) 



27rfc 2 



"»,(£, +nii) 



(D6) 



The spin factor S SlS2 - s ' s > of the interaction can be writ- 
ten as 
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s s 1 s 2 - s [s' 2 = [u(k, si)7o(l)u(fe', s' x )] [v(-k, s 2 )7o(2)u(-fe', 4)] ~ K fc , Si)ji(i)u(k' , s[)} [v(-k, s 2 )ji(2)v(-k' , 4)] 

+ [^J(fc)B 1 (fc / )«T(l)a- fc /(l) + Br(AMi(fc / )a- fc (l)«T(l)] 

■[i? 2 *(fc)A 2 (fc'K(2)^(2)+^(fc)i?2(fc')^(2)^(2)]}|x(s' 1 )x(4))- (D7) 



1 / 



47T 



(D8) 



The kernel Ugi. s i V {k] k') can be rewritten as 

Ui tS , l ,(k;k') = (^j s iM(n k ,s)\U[k,k';<T(l),<T(2)]\^ Js , VM (^,s')), 
where the interaction operator in momentum and spin space is 

+ [A*(fc) J B 1 (fc')T(l)^'(l) + S 1 *(fc)A 1 (fc'K(l)T(l)] 

• [B|(fc)A 2 (fc')^fc(2)«T(2) + A* 2 (k)B 2 (k')a(2)a k ,(2)} }. (D9) 

I 

Since a^ and ov are pseudo scalar, k, k', a^ ov, and and triplet. 

<r(l) -^(2) arc scalar, the above interaction kernel oper- Jf ^ = ^ and fche tengor potcntials are neglectedj , 

ator £/[«, k ; <t(1J,<t(/JJ is scalar. an( j g are conse rved and the interaction kernel becomes 

From the last expression of the kernel U J sl . s , h , (k; k'), we diagonal in I and s representation, 
could see that the first term contributes to different kinds 
of central potentials and relativistic spin-orbit coupling 

potentials, the second term contributes to the tensor po- U^ l . s , l ,{k; k') = U^. sl (k; k')5u>5 ss > — Uj s i(k; k')5u>5 ss i. 

tentials changing I by Al — ±2 and mixing spin singlet (D10) 
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